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Microbial  fuel  cells  (MFCs)  offer  great  promise  for  simultaneous  treatment  of  wastewater  and  energy 
recovery.  While  past  research  has  been  based  extensively  on  experimental  studies,  modelling  and  sim¬ 
ulation  remains  scarce.  A  typical  MFC  shares  many  similarities  with  chemical  fuel  cells  such  as  direct 
ascorbic  acid  fuel  cells  and  direct  methanol  fuel  cells.  Therefore,  an  attempt  is  made  to  develop  a  MFC 
model  similar  to  that  for  chemical  fuel  cells.  By  integrating  biochemical  reactions,  Butler-Volmer  expres¬ 
sions  and  mass/charge  balances,  a  MFC  model  based  on  a  two-chamber  configuration  is  developed  that 
simulates  both  steady  and  dynamic  behaviour  of  a  MFC,  including  voltage,  power  density,  fuel  concen¬ 
tration,  and  the  influence  of  various  parameters  on  power  generation.  Results  show  that  the  cathodic 
reaction  is  the  most  significant  limiting  factor  of  MFC  performance.  Periodic  changes  in  the  flow  rate  of 
fuel  result  in  a  boost  of  power  output;  this  offers  further  insight  into  MFC  behaviour.  In  addition  to  a  MFC 
fuelled  by  acetate,  the  present  method  is  also  successfully  extended  to  using  artificial  wastewater  (solu¬ 
tion  of  glucose  and  glutamic  acid)  as  fuel.  Since  the  proposed  modelling  method  is  easy  to  implement,  it 
can  serve  as  a  framework  for  modelling  other  types  of  MFC  and  thereby  will  facilitate  the  development 
and  scale-up  of  more  efficient  MFCs. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  microbial  fuel  cell  (MFC)  is  a  device  capable  of  directly  trans¬ 
forming  chemical  energy  into  electrical  energy  via  electrochemical 
reactions  involving  biochemical  pathways.  Unlike  conventional  fuel 
cells,  which  usually  use  precious  metals  as  catalysts,  MFCs  use  living 
organisms  at  mild  conditions,  namely,  room  temperature,  atmo¬ 
spheric  pressure,  and  neutral  pH  [1].  MFCs  have  operational  and 
functional  advantages  over  the  technologies  currently  used  for  gen¬ 
erating  energy  from  organic  matter  [2].  One  of  the  most  active 
areas  of  MFC  research  in  the  past  few  years  has  been  the  produc¬ 
tion  of  power  from  organic  waste  [3].  While  full  scale  and  highly 
effective  MFCs  are  not  yet  available,  the  technology  holds  great 
promise.  The  growing  pressure  for  sustainable  development  and 
the  call  for  renewable  energy  further  motivate  the  development  of 
this  technology.  Over  the  past  years,  considerable  effort  has  been 
made  to  study  the  biological  aspects  such  as  the  molecular  analy¬ 
sis  of  the  microbial  community  and  electron  transfer  to  electrodes 
in  microorganisms  [4,5],  as  well  as  developing  various  designs  and 
configurations  of  MFCs  [6,7].  Recent  studies  on  the  conversion  of 
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organic  wastes  to  electricity  have  focused  on  the  engineering  chal¬ 
lenges  [3]. 

A  number  of  limitations  have  hindered  the  wide  implementa¬ 
tion  of  MFCs.  In  particular,  the  power  density  of  a  MFC  is  several 
orders  of  magnitude  lower  than  that  of  chemical  fuel  cells,  and 
the  technology  is  still  only  in  the  laboratory  phase  [8,9].  The 
extractable  power  from  a  particular  fuel  cell  is  restricted  by  a 
range  of  different  parameters  [2]  that  include  the  amount  of  bac¬ 
terial  cells,  mixing  and  mass  transfer  phenomena  in  the  reactors, 
bacterial  kinetics,  cathodic  reactions,  and  the  efficiency  of  the 
proton-exchange  membrane.  Therefore,  construction  and  analy¬ 
sis  of  MFCs  involve  multidisciplinary  knowledge  of  microbiology, 
electrochemistry,  materials  science,  and  engineering.  Various  intra¬ 
system  phenomena  and  the  effects  of  operating  factors  are  not  well 
understood.  Investigating  each  parameter  independently  and/or  all 
possible  combinations  of  parameters  using  laboratory  experiments 
is  both  costly  and  time  consuming.  Numerical  modelling  of  a  fuel 
cell  system  is  a  valuable  tool  for  investigating  system  parameters 
with  reduced  time  and  money,  as  models  can  be  easily  modified 
to  simulate  various  configurations  and  operating  conditions.  The 
development  of  a  successful  commercial  MFC  requires  a  ‘chemical 
engineering’  approach  [10].  While  modelling  and  simulation  has 
been  widely  used  to  develop  various  chemical  fuel  cell  systems, 
modelling  and  other  forms  of  quantifications  are  often  overlooked 
in  MFC  development  [11  ].  Since  the  first  MFC  model  was  proposed 


0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2009.06.101 


80 


Y.  Zeng  et  al.  /  Journal  of  Power  Sources  195  (2010)  79-89 


[12]  to  simulate  a  suspended  cell  with  an  added  redox  mediator, 
very  little  progress  has  been  made  in  this  field.  Only  very  recently 
have  computational  models  for  biofilm-based  microbial  fuel  cells 
been  reported  [  13,14].  Since  the  models  are  comprehensive  because 
of  detailed  descriptions  of  multi-species  for  the  biofilm  anode,  mod¬ 
elling  methods  may  not  be  readily  implemented  by  the  majority  of 
the  MFC  community. 

A  typical  two-chamber  microbial  fuel  cell  shares  some  similar¬ 
ities  with  chemical  fuel  cells  such  as  direct  methanol  fuel  cells 
(DMFCs)  [15]  and  direct  ascorbic  acid  fuel  cells  (DAAFCs)  [16]. 
In  both  cases,  anodic  and  cathodic  chambers  are  separated  by  a 
cation-exchange  membrane,  fuels  are  oxidized  in  the  anodic  cham¬ 
ber  and  release  protons  that  pass  through  the  membrane  to  reach 
the  cathodic  chamber.  Water  saturated  with  oxygen  is  fed  into  the 
cathodic  chamber,  where  the  transported  protons  combine  with 
dissolved  oxygen.  A  MFC  system  is,  however,  remarkably  differ¬ 
ent  from  chemical  fuel  cells  in  terms  of  the  anodic  reaction,  fuel 
concentration,  ion  strength,  pH,  operating  temperature,  and  purity 
of  fuels.  In  particular,  the  fuel  oxidation  in  a  chemical  fuel  cell  is 
often  catalyzed  by  noble  metals  under  high  temperature,  whereas 
a  wide  range  of  organic  compounds  can  be  catalyzed  in  a  MFC  by 
micro-organisms  at  room  temperature.  In  addition,  a  variety  of  inor¬ 
ganic  chemicals  are  employed  to  support  the  microbial  metabolism 
in  a  MFC  [5,17,18].  As  a  result,  much  more  complicated  bio¬ 
electrochemical  reactions  often  take  place  in  the  anodic  chamber  of 
a  MFC.  Despite  of  the  differences  between  chemical  fuel  cells  and 
MFCs,  we  have  been  motivated  by  the  successful  development  of 
mathematical  models  of  a  DMFC  [15]  and  a  DAAFC  [16]  that  config¬ 
urations  similar  to  that  of  a  typical  two-chamber  MFC.  In  the  present 
study,  these  methodologies  have  been  adopted  to  develop  the  MFC 
model  by  taking  into  account  the  bio-electrochemical  reactions. 

In  the  anodic  compartment  of  a  MFC,  biological  and  electro¬ 
chemical  reactions  take  place  via  various  mechanisms  of  electron 
transfer  such  as  biofilm  [14]  and  planktonic  [19].  In  this  work,  the 
model  formulation  is  simplified  with  lump  parameters  of  the  fuel 
oxidation.  Since  the  anodic  chamber  of  a  MFC  operates  at  anaer¬ 
obic  or  anoxic  conditions,  the  biochemical  kinetics  are  simulated 
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Fig.  1.  Schematic  diagram  of  acetate  MFC  experiment. 


in  a  similar  fashion  to  that  of  a  Anaerobic  Digestion  Model  No.  1 
(ADM1)  [20].  As  in  the  case  of  DAAFC,  the  rate  of  reaction  is  con¬ 
trolled  by  the  electrical  potential  in  an  electrochemical  cell  and 
therefore  the  Butler-Volmer  expression  is  incorporated.  In  compar¬ 
ison  with  chemical  fuel  cells,  which  use  only  pure  fuel,  MFCs  often 
employ  a  variety  of  inorganic  chemicals  to  support  the  microbial 
metabolism  in  addition  to  the  fuel,  which  causes  cations  other  than 
protons  to  dominate  the  mass  transport  via  the  proton-exchange 
membrane  [5,17,18].  This  issue  is  also  addressed  in  the  MFC  mod¬ 
elling.  As  a  first  attempt,  a  mediator-less  two-compartment  MFC 
with  a  cation-exchange  membrane  was  studied,  and  acetate  was 
fed  as  the  fuel  since  it  is  the  most  abundant  fatty  acid  in  anaerobic 
ecosystems  and  is  used  as  an  electron  donor  by  anaerobic  respira¬ 
tory  bacteria.  The  modelling  method  was  extended  to  an  artificial 
wastewater  MFC  fed  with  solution  of  glucose  and  glutamic  acid. 

The  model  is  developed  by  integrating  the  bio-electrochemical 
kinetics  and  mass  and  charge  balances  within  the  MFC.  Various 
parameters  that  significantly  affect  the  MFC  performance  are  inves¬ 
tigated  in  the  simulation,  and  the  evolution  in  time  of  the  MFC 
voltage,  power  density  and  fuel  concentration  are  studied  through 
dynamic  simulation  with  stepped  and  periodic  changes  in  the 
anodic  feed  flowrate.  The  developed  method  of  MFC  modelling  is 
relatively  simple  and  easy  to  implement,  which  may  enhance  its 
practical  use  in  the  design  and  operation  of  MFCs. 

2.  Methods  of  approach 

2.1.  Experiments 

A  schematic  diagram  of  the  experimental  set-up  is  shown  in 
Fig.  1.  The  working  principle  of  the  microbial  fuel  cell  can  be  found 
elsewhere  [21,22].  In  the  most  general  sense,  microbial  fuel  cells 
function  by  oxidizing  an  electron  donor  with  electron  transfer  at 
the  anode  under  anoxic  conditions.  Electrons  donated  by  the  anode 
pass  through  a  resistor  or  other  type  of  electrical  device  to  the  cath¬ 
ode,  which  is  submerged  in  aerobic  water.  In  this  work,  the  MFC 
device  was  made  of  transparent  polyacrylic  plastic.  Each  MFC  con¬ 
sisted  of  anode  and  cathode  compartments  of  equal  volume  (5.5  ml) 
and  dimensions  (5.5  x  1.0  x  1.0  cm).  Each  compartment  contained 
two  pieces  (4.5  x  1.0  x  0.5  cm)  of  graphite  felt  (GF  series,  Electro¬ 
synthesis,  Amherst,  NY)  as  electrodes,  but  the  graphite  felt  used 
for  the  cathode  electrode  was  coated  with  0.3  mg  cm-2  of  platinum 
powder  [23  ].  The  electrodes  were  connected  via  a  platinum  wire  to  a 
resistance  box  and  a  multimeter  (Keithley  Instruments,  Cleveland, 
OH,  USA).  The  anode  and  cathode  chambers  were  separated  by  a 
cation-exchange  membrane  (Nafion®,  Dupont  Co.,  USA).  Injection 
ports  were  installed  in  each  chamber  of  the  fuel  cell.  The  anode 
chamber  was  kept  anoxic  by  purging  nitrogen  gas.  Throughout 
the  study,  the  fuel  containing  1.56  mM  acetate  and  buffer  solution 
(4mM  sodium  acetate  and  inorganic  salts)  was  continuously  fed 
to  the  anode  at  the  rate  of  0.375  ml  min-1,  and  air-saturated  tap 
water  was  fed  into  the  cathode  chamber  as  oxidant  at  the  rate  of 
18.5  ml  min-1.  Peristaltic  pumps  (505S,  Watson-Marlow,  Campel, 
UK)  were  used  to  feed  the  liquids  that  were  pre-warmed  to  30  °C. 
All  experiments  were  conducted  using  three  parallel  microbial  fuel 
cells  installed  in  a  temperature-controlled  chamber  at  30  °C  and 
operated  over  one  year.  The  external  resistance  between  the  anode 
and  the  cathode  was  systematically  varied  by  means  of  a  resis¬ 
tance  box  in  the  range  of  10-10,000  £2.  The  potential  between  the 
anode  and  the  cathode  was  measured  with  the  multimeter  and 
recorded  on  a  personal  computer  through  a  data-acquisition  sys¬ 
tem  (Testpoint,  Capital  Equipment,  Richmond,  VA,  USA).  The  data 
was  taken  only  when  steady-state  conditions  had  been  established 
after  changing  the  external  resistance  and  took  up  to  30  min.  The 
average  or  typical  readings  of  the  three  fuel  cells  were  used  in  model 
development  and  verification. 
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2.2.  Model  development 


Experiments  suggest  that  acetate  is  oxidized  in  the  anode  by  the 
reaction  of  an  eight-electron  transfer  as  described  in  Eq.  (R1 ).  Since 
the  anode  chamber  operates  under  anaerobic  or  anoxic  conditions, 
the  acetate  degradation/oxidation  can  be  modelled  in  a  similar 
manner  to  ADM1,  where  a  Monod-type  equation  is  adopted.  In  addi¬ 
tion,  the  operating  conditions  remain  nearly  unchanged  except  for 
variations  in  the  external  resistances.  Various  microbial  consortia 
are  thus  lumped  to  a  single  quantity  (biomass).  Since  (Rl)  is  a  bio¬ 
electrochemical  reaction  controlled  by  the  electrical  potential  in 
an  electrochemical  cell,  the  Butler-Volmer  expression  is  incorpo¬ 
rated.  Furthermore,  since  the  reverse  reaction  is  insignificant,  only 
the  forward  reaction  is  used  to  formulate  the  reaction  rate  (Eq.  (1 )). 
It  is  noted  that  Eq.  (1)  is  similar  to  the  expression  in  the  literature 
[14]. 

(CH20)2  +  2H20  ->■  2C02  +  8H+  +  8e“  (Rl ) 

where:  CAc  and  X  are  the  concentrations  of  acetate  and  biomass 
in  the  anode  compartment,  respectively;  pa  is  the  anodic  overpo¬ 
tential;  is  the  rate  constant  of  the  anode  reaction  at  standard 
conditions  (maximum  specific  growth  rate);  KAc  is  the  half  velocity 
rate  constant  for  acetate;  a  is  the  charge  transfer  coefficient  of  the 
anodic  reaction,  F  is  the  Faraday  constant,  R  is  the  gas  constant,  T  is 
the  cell  operating  temperature. 

Since  MFCs  employ  a  variety  of  inorganic  chemicals  to  support 
the  microbial  metabolism,  cations  such  as  I<+,  Na+,  NH4+,  Mg2+, 
and  Ca2+  are  dissociated  and  their  concentrations  are  typically  105 
times  higher  than  that  of  protons  at  neutral  pH.  Thus  the  number  of 
cations  transported  from  the  anode  to  the  cathode  compartments 
other  than  protons  are  the  same  as  the  number  of  electrons  trans¬ 
ferred  through  the  circuit  [5,17,18].  This  indicates  that  virtually  no 
protons  are  transported  in  the  MFC,  and  electro-neutrality  is  sus¬ 
tained  mainly  by  the  transport  of  cations  instead  of  protons.  An 
analysis  [17]  of  a  Nation  117  membrane  in  a  two-compartment  MFC 
shows  that  K+  and  Na+  occupied  about  74%  of  the  sulfonate  residues. 
In  order  to  maintain  Nation  conductivity  in  the  long  run,  the  Nation 
membrane  was  boiled  in  0.1  M  HC1  before  use  to  replace  any  cations 
attached  to  the  sulfonate  residues  with  protons.  It  is  believed  that 
the  cation  occupation  is  not  permanent,  but  transient.  Under  these 
conditions,  it  is  still  possible  for  the  cathode  reaction  to  consume 
protons  available  from  the  dissociation  of  water. 

For  the  sake  of  simplicity,  it  is  assumed  that  only  univalent 
cations  M+  transport  through  the  membrane  and  the  M+  ions  do 
not  involve  in  the  reaction  at  cathode.  Consequently,  the  reduction 
of  dissolved  oxygen  in  the  cathode  is  suggested  as: 

02  +  4e“  +  2H20  ->  40H-  (R2) 


It  was  found  that  dissolved  oxygen  exhibits  Monod-type 
behaviour  [5]  and  our  preliminary  study  shows  that  the  reverse 
reaction  of  oxygen  reduction  is  negligible.  Therefore,  the  rate  of 
reaction  in  the  cathodic  chamber  is  formulated  as  Eq.  (2).  Again, 
the  Butler-Volmer  expression  is  incorporated  to  describe  the  elec¬ 
trochemical  reaction. 


r2 


-4 


Co2 

Kq2  +  C02 


exp 


(2) 


where:  C02  is  the  concentration  of  dissolved  oxygen  in  the  cathode 
compartment;  rjc  is  the  overpotential  at  the  cathode;  Kq2  is  the  half¬ 
velocity  rate  constant  for  dissolved  oxygen;  is  the  rate  constant 
of  the  cathode  reaction  under  standard  conditions;  j3  is  the  charge- 
transfer  coefficient  of  the  cathodic  reaction.  Water  concentration  is 
assumed  constant  (excess  component  in  liquid  mixture). 


Similar  to  DAAFC  modelling  [16],  it  is  assumed  that  both  the 
anode  and  cathode  compartments  can  be  treated  as  a  continuously 
stirred  tank  reactor  (CSTR).  Phase  mixture,  i.e.,  all  mass-transport 
processes,  is  assumed  to  be  so  fast  compared  with  the  biochemical 
and  redox  reactions,  such  that  the  concentrations  of  all  reactants  in 
the  bulk  solution  can  be  regarded  to  be  equal  to  those  on  the  surface 
of  electrodes.  In  addition,  carbon  dioxide  and  acetate  are  assumed 
not  to  diffuse  into  the  membrane,  and  the  gas-phase  formation  by 
release  of  carbon  dioxide  bubbles  is  not  taken  into  account.  Conse¬ 
quently,  the  mass  balances  of  the  four  components  in  the  anode 
compartment,  namely,  acetate,  dissolved  C02,  hydrogen  ion  and 
biomass,  are  expressed  by  Eqs.  (3)-(6),  respectively: 


d<ff  =  Qa(C‘Anc-CAc)-Amri 

(3) 

=  Qa(C'n02  -  CCo2 )  +  2 Amr, 

(4) 

^L=Qa(ClHn-CH)  +  8Amn 

(5) 

=Qa(X"yX)  +AmYacri  VaKdecX 

(6) 

In  the  above  equations,  the  subscripts  ‘a’  and  ‘in’  denote  the  anode 
and  the  feed  flow,  respectively.  V,  Q  and  Am  are  the  volume  of 
the  compartment,  the  feed  flow  rate,  and  the  cross-section  area 
of  membrane,  respectively.  In  Eq.  (6),  fx  represents  the  reciprocal 
of  the  wash-out  fraction,  Yac  the  bacterial  yield,  and  Kdec  the  decay 
constant  for  acetate  utilisers. 

In  the  cathode  compartment,  the  mass  balances  of  dissolved  02, 
hydroxyl,  and  cation  M+  are  expressed  by  Eqs.  (7)-(9),  respectively: 


Vc^=Qc(C",2-c°2)  +  ryim 

(7) 

=  Qc(Cq^j  -  C0H)  -  4r2A,„ 

(8) 

vc^fL=QActi-CM)  +  NMAm 

O) 

The  subscript  ‘c’  denotes  the  cathode.  In  Eq.  (9),  NM  is  the  flux  of 
M+  ions  transported  from  the  anode  to  cathode  compartment  via 
the  membrane.  It  is  noted  that  the  following  relationship  is  held  for 
the  cell  current  density  and  the  flux  of  ions  via  the  membrane  [24]: 

icell=Fj2^i  (10) 


where:  z*  is  the  charge  number  of  the  ith  species;  N{  is  the  superfi¬ 
cial  flux  of  the  ith  species,  ice //  denotes  the  cell  current  density.  As 
mentioned  previously,  only  M+  ions  are  assumed  to  be  transported 
through  the  membrane,  thus  the  flux  of  M+  ion  (mol  m-2  h-1 )  can 
be  calculated  as  follows,  where  the  coefficient  3600  is  the  factor  of 
unit  conversion. 

JVm=  3600i££|,  (11) 

The  charge  balances  at  the  anode  and  cathode  are  given  by  Eqs. 
(12)  and  (13),  respectively,  where  Ca  and  Cc  are  the  capacitances  of 
the  anode  and  cathode,  respectively. 

Cad^-=3600icell-8Fr]  (12) 

Cc^  =  -3600iceu-4Fr2  (13) 

It  is  assumed  that  the  ohmic  drops  in  the  current-collectors 
and  electric  connections  are  negligible,  and  the  cell  resistance  is 
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Table  1 

Parameters  of  experiment  set-up  and  operation,  and  constants  of  acetate  MFC  model. 


Symbol 

Description 

Unit 

Value 

F 

Faraday’s  constant 

Coulombs  mol-1 

96485.4 

R 

Gas  constant 

J  mol-1  K-1 

8.3144 

T 

Temperature 

K 

303 

km 

Electrical  conductivity  of  membrane 

Ohm-1  m1 

17 

dm 

Thickness  of  membrane 

m 

1.778  x  10“4 

kaQ 

Electrical  conductivity  of  the  aqueous  solution 

Ohm-1  m1 

5 

(JceH 

Distance  between  anode  and  cathode  in  the  cell 

m 

2.2  x  10-2 

Ca 

Capacitance  of  anode 

Fm-2 

4x  102 

Cc 

Capacitance  of  cathode 

Fm-2 

5  x  102 

Va 

Volume  of  anode  compartment 

m3 

5.5  x  10“5 

Vc 

Volume  of  cathode  compartment 

m3 

5.5  x  10“5 

Am 

Area  of  membrane 

m2 

5  x  10-4 

Yac 

Bacterial  yield 

Dimensionless 

0.05 

Kdec 

Decay  constant  for  acetate  utilisers 

h-1 

8.33  x  10-4 

fx 

Reciprocal  of  wash-out  fraction 

Dimensionless 

10 

On 

Flow  rate  of  fuel  feed  to  anode 

m3  h_1 

2.25  x  10“5 

Qc 

Flow  rate  feeding  to  cathode  compartment 

m3  h_1 

1.11  x  10“3 

fin 
^ Ac 

Concentration  of  acetate  in  the  influent  of  anode  compartment 

mol  ITT3 

1.56 

fin 

co2 

Concentration  of  C02  in  the  influent  of  anode  compartment 

mol  m3 

0 

)Cn 

Concentration  of  bacteria  in  the  influent  of  anode  compartment 

mol  m3 

0 

fin 
^ H 

Concentration  of  H+  in  the  influent  of  anode  compartment 

mol  m3 

0 

fin 

o2 

Concentration  of  dissolved  02  in  the  influent  of  cathode  compartment 

mol  m3 

0.3125 

fin 

U M 

Concentration  of  M+  in  the  influent  of  cathode  compartment 

mol  m3 

0 

fin 

COH 

Concentration  of  OH  in  the  influent  of  cathode  compartment 

mol  m3 

0 

U° 

Cell  open  circuit  potential 

volt 

0.77 

solely  due  to  the  resistances  of  the  membrane  and  the  solution. 
Consequently,  the  cell  voltage  Uceu  is  calculated  as: 

U«fl  =  U°-i?a  +  »;c-(^  +  ^)fceB  (14) 

where:  U°  is  the  open-circuit  voltage;  dm  and  dce//  are  the  thickness 
of  the  membrane  and  the  distance  of  the  electrodes,  respectively; 
km  and  kaq  are  the  conductivities  of  the  membrane  and  the  solution, 
respectively. 

2.3.  Parameter  estimation 

For  the  model  system  of  Eqs.  (1)-(14),  the  operating  conditions, 
the  parameters  of  the  experimental  set-up  and  the  constants  are 
listed  in  Table  1.  Among  them,  the  electrical  conductivity  of  the 
aqueous  solution  ( kaq )  is  estimated,  but  it  is  not  sensitive  to  simu¬ 
lation  results.  The  parameters  of  the  capacitances  of  the  anode  and 
cathode  are  not  required  for  steady-state  simulation  and  parameter 
estimation,  but  they  are  required  to  study  the  dynamic  simulation. 
Similarly,  although  the  two  parameters  are  estimated  from  another 
study,  they  are  not  sensitive  to  the  present  results.  Both  bacterial 
yield  and  decay  constant  for  acetate  utilisers  are  taken  from  [25]. 

The  six  model  parameters  listed  in  Table  2  cannot  be  directly 
determined  from  experiments  and  thus  are  estimated  by  a  mathe¬ 
matical  method  of  best  fitting  of  the  experimental  data.  Specifically, 


the  model  parameters  are  estimated  by  minimizing  the  absolute 
differences  between  the  measured  and  simulated  cell  voltages, 
namely,  minimizing  the  quantity  of  Fobj  defined  in  Eq.  (15). 

N 

(15) 

1  =  1 

where  N  is  the  number  of  data  points,  U^xp  is  the  experimental 
measurements,  and  the  simulated  cell  voltage  U Fal  is  obtained  from 
Eq.  ( 14)  by  solving  the  model  system  of  Eqs.  ( 1 )-( 13 )  at  the  steady- 
state,  namely,  the  algebraic  equations  obtained  by  setting  all  the 
derivatives  to  zeros. 

3.  Results  and  discussion 

As  mentioned  early,  three  identical  acetate  MFC  devices  were 
used  to  conduct  parallel  experiments.  During  stable  operation,  the 
external  resistance  (Rex)  was  varied  and  then  fixed  at  17  points  in 
the  range  of  10-10,000  £2.  For  each  resistance,  the  MFC  was  allowed 
to  establish  its  steady-state,  which  took  up  to  30  min.  The  resulting 
cell  voltages  (Uceu)  were  recorded  and  the  cell  currents  (JceH)  were 
calculated  with  Ice[[  =  Uceu/Rex.  Each  data  point  was  obtained  by  tak¬ 
ing  the  average  readings  of  the  three  devices,  and  17  data  points  in 
total  were  obtained.  Among  them,  10  data  points  were  selected  for 


Table  2 

Estimated  parameters  of  acetate  MFC. 


Symbol 

Description 

Unit 

Resultant  values 

k] 

Forward  rate  constant  of  anode  reaction  at  standard  condition  (maximum  specific  growth  rate) 

mol  irr2  h-1 

0.207 

k°2 

Forward  rate  constant  of  cathode  reaction  at  standard  condition 

m12  mol-4  h_1 

3.288  x  10-5 

Kac 

Half  velocity  rate  constant  for  acetate 

mol  ITT3 

0.592 

Ko2 

Half  velocity  rate  constant  for  dissolved  oxygen 

mol  ITT3 

0.004 

a 

Charge  transfer  coefficient  of  anode 

Dimensionless 

0.051 

V 

Charge  transfer  coefficient  of  cathode 

Dimensionless 

0.663 
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Fig.  2.  Results  from  fitting  experimental  data  (a)  and  model  validation  (b)  for  acetate  MFC. 


Table  3 

Percentage  changes  in  power  density  of  acetate  MFC  when  each  parameter  value  is  varied  one  at  a  time. 


Ratio 

*? 

k° 

I<Ac 

I< o2 

Of 

P 

Qa 

Qc 

cin 

U Ac 

q?2 

£ pen 

dm 

kaq 

0.8 

-43.00 

-6.55 

18.37 

0.08 

10.23 

63.01 

-10.57 

-0.01 

-78.42 

-0.12 

1.86 

0.01 

-2.32 

1.2 

35.13 

5.35 

-16.71 

-0.08 

-6.82 

-144.71 

-17.47 

0.01 

31.09 

0.08 

-1.86 

-0.01 

1.55 

parameter  estimation  by  best  fitting  the  model,  and  the  remaining 
7  data  were  used  for  model  validation. 


3.1  Resultant  parameters  and  sensitivity  analysis 

As  shown  in  Fig.  2a,  good  agreement  is  observed  between  the 
calculated  and  the  measured  cell  voltages  for  a  fitting  set  with  the 
use  of  the  6  resultant  model  parameters.  The  average  absolute  dif¬ 
ference  is  0.01  V.  The  model  is  further  examined  by  a  validation 
set.  Again,  good  agreement  is  found  as  shown  in  Fig.  2b  where  the 
average  absolute  difference  is  0.01  V. 

A  preliminary  sensitivity  analysis  on  the  operating,  design  and 
model  parameters  was  investigated  by  varying  each  parameter,  one 
at  a  time,  while  leaving  the  other  parameters  unchanged.  Specif¬ 
ically,  each  parameter  is  varied  by  multiplying  the  ‘ratio’  to  its 
original  value  at  the  base  case— the  experimental  condition.  The 
resultant  average  percentage  changes  in  the  power  output  over  the 
range  of  cell  current  density  (1-11.5  Am-2)  are  listed  in  Table  3. 
As  shown  in  the  Table,  when  the  parameter  is  reduced  by  a  factor 
of  0.8,  the  acetate  feed  concentration  leads  to  the  largest  decrease 
in  power  density  over  the  entire  range  of  cell  current  density,  i.e., 
there  is  an  absolute  change  of  78%.  The  next  largest  change  is  the 
electron  transfer  coefficient  of  the  cathode  (ft),  which  leads  to  a 
boost  of  63%.  On  the  other  hand,  when  the  parameter  is  increased 
by  a  factor  of  1.2,  the  cathode  electron  transfer  coefficient  results  in 
the  largest  change,  i.e.,  a  decrease  of  114%.  The  second  and  the  third 
largest  changes  are  caused  by  the  rate  constant  of  the  anode  reac¬ 
tion  (35%)  and  the  acetate  feed  concentration  (31%).  In  both  cases, 
the  changes  in  the  electrical  conductivity  of  the  aqueous  solution 
kaq  is  very  small,  namely,  about  2%,  which  implies  that  it  is  insignif¬ 
icant  to  the  accuracy  of  model  prediction.  It  is  also  noted  that  the 
impacts  of  both  cathode  feed  flowrate  and  dissolved  oxygen  con¬ 
centration  are  very  insignificant,  and  membrane  thickness  is  the 
least  sensitive. 

The  sensitivity  of  the  6  model  parameters  is  further  studied  by 
the  local  relative  sensitivity  analysis  method  [26],  to  evaluate  the 
ratio  of  changes  in  the  computed  power  density  to  the  changes  in 
the  parameters.  The  following  equation  is  used  for  the  6  model 
parameters. 


P(t,  Xj  +  Sxj)  -  P(t ,  Xj) 

8xj 


P(t,XjY 


where:  7)  is  time-dependent  sensitivity  of  the  jth  parameters;  Xj 
is  the  value  of  jth  parameter;  8xj  is  the  change  in  Xj\  and  P  is  the 
power  density.  In  the  present  study,  8xj  =  0.0\Xj.  A  dynamic  sim¬ 
ulation  of  step  change  in  acetate  feed  flowrate  from  2.2  x  10-5 
to  1  x  10-5  (m3h-1)  was  conducted  to  examine  the  sensitivities 
of  the  6  parameters.  The  results  are  shown  in  Fig.  3,  and  the 
order  of  parameter  sensitivity  (from  the  highest  to  the  lowest)  is 
/3>k(l>a>KAc>l<°>Ko2.  Apparently,  the  electron  transfer  coeffi¬ 
cient  is  the  most  sensitive  parameter,  whereas  I<q2  is  the  least 
significant  and  remains  almost  unchanged. 

The  parameter  sensitivity  analysis  has  shown  that  the  cathodic 
reaction  may  be  the  most  significant  factor  limiting  the  perfor¬ 
mance  of  MFCs.  This  result  agrees  with  the  findings  of  Zhao  et  al. 
[  18],  who  reported  that  the  cathodic  reaction  often  limits  the  perfor¬ 
mance  of  chemical  fuel  cells,  such  as  polymer  electrolyte  membrane 
fuel  cells  and  solid  oxide  fuels,  and  MFCs  also  share  this  problem.  To 
improve  the  cathodic  reaction,  one  of  the  most  important  issues  is  to 
develop  efficient  cathode  materials  which  is,  however,  beyond  the 
scope  of  this  study  and  should  be  the  focus  of  further  investigation. 

3.2.  Steady-state  simulation 

With  the  parameters  listed  in  Table  2,  the  component  concentra¬ 
tions  and  reaction  rates  are  evaluated  with  respect  to  the  changes 


X 


J  =  1 . 6 


(16) 


Fig.  3.  Resultant  sensitivities  of  6  model  parameters  of  acetate  MFC  model. 
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Fig.  4.  Results  of  steady-state  simulation  of  acetate  MFC.  (a)  Reaction  rates  of  acetate  oxidation  n  and  oxygen  reduction  r2;  (b)  concentrations  of  acetate,  dissolved  CO2, 
biomass  (X),  dissolved  oxygen  (02)  and  M+;  (c)  pH  values  in  anode  and  cathode  compartment;  (d)  overpotentials  of  anode  and  cathode. 


in  cell  current  density.  As  shown  in  Fig.  4a,  the  reaction  rate  of 
acetate  oxidation  in  the  anode  chamber  (ri)  and  the  rate  of  oxy¬ 
gen  reduction  (-r2)  in  the  cathode  chamber  are  proportional  to  the 
cell  current  density.  This  behaviour  can  be  explained  by  Eqs.  (12) 
and  (13)  at  the  steady-state  that  are  reduced  to  r\  =450 iceulF  and 
r2  =  -900 iceulF,  respectively. 

The  resultant  concentration  of  5  components,  i.e.,  acetate,  C02, 
biomass  (X)  in  the  anodic  chamber,  and  dissolved  oxygen  and  M+ 
in  the  cathodic  chamber,  are  plotted  in  Fig.  4b.  Obviously,  C02  and 
biomass  linearly  increase  while  acetate  decreases  with  respect  to 
iceii.  Indeed,  a  higher  current  density  leads  to  a  faster  reaction  which 
consumes  more  fuel  (acetate)  and  generates  more  C02.  Similarly,  a 
higher  current  density  results  in  faster  decrease  in  the  dissolved 
oxygen  in  the  cathode  chamber  although  it  is  not  pronounced.  It 
follows  from  Eq.  (9)  that  M+  concentration  is  proportional  to  the 
M+  flux,  which  is  a  linear  function  of  cell  current  density  (Eq.  (11 )). 
Therefore,  M+  concentration  linearly  increases  with  current  den¬ 
sity. 


Fig.  5.  Simulation  results  of  variation  in  (a)  acetate  feed  flowrate  (Qa);  (b)  acetate 
of  0.8  and  1.2,  respectively. 


It  can  be  seen  in  Fig.  4c  that  the  pH  in  the  anode  compart¬ 
ment  monotonically  decreases  as  the  cell  current  density  increases, 
which  indicates  the  anode  compartment  is  acidified  and  the  anode 
potential  is  raised  (Fig.  4d).  This  result  is  in  agreement  with  the 
experiment  observations  by  other  research  groups  [8,21].  By  con¬ 
trast,  the  pH  in  the  cathodic  compartment  increases  monotonously, 
which  also  agrees  with  the  phenomena  reported  in  [5]. 

As  revealed  from  the  parameter  sensitivity  analysis,  the  acetate 
flowrate  (Qa)  and  concentration  (C™)  in  the  feed  flow  demon¬ 
strate  the  highest  sensitivities  among  the  operational  and  design 
parameters.  The  effects  of  the  two  parameters  have  therefore  been 
further  examined.  Specifically,  the  two  parameters  have  been  var¬ 
ied  by  multiplying  their  original  values  by  a  factor  of  0.8  and  1.2, 
respectively,  while  the  other  parameters  remain  unchanged.  The 
simulation  results  are  compared  with  that  of  the  original  experi¬ 
mental  set-up  (the  Base  Case).  As  shown  in  Fig.  5a,  a  decrease  in 
acetate  feed  flowrate  (Qa)  leads  to  a  boost  in  power  output  until  the 
cell  current  density  is  higher  than  8.5  Am-2.  Beyond  this  point,  an 


concentration  (C).  Qa  and  C  are  varied  by  multiplying  their  original  values  by  a  factor 
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Fig.  6.  Results  of  dynamical  simulation  of  step  change  in  acetate  feed  flowrate  from  2.25  x  10-5  to  1  x  10-5  (m3  h-1):  (a)  acetate  concentration;  (b)  anode  overpotential;  (c) 
cell  voltage;  (d)  power  density. 


extremely  sharp  decrease  in  power  output  is  observed.  An  increase 
in  Qa  demonstrates  opposite  behaviour.  The  result  implies  that  a 
lower  fuel  feed  flowrate  may  increase  the  power  output  but  at  the 
expense  of  a  reduction  in  attainable  current  density.  In  order  to 
maintain  the  same  range  of  current  density,  the  fuel  concentration 
has  to  be  increased.  The  effect  of  acetate  feed  concentration  C™c  is 
shown  in  Fig.  5b.  An  opposite  effect  to  that  of  the  feed  flowrate 
is  observed,  i.e.,  cell  current  density  is  boosted  by  an  increase 
in  the  feed  concentration  and  the  attainable  current  density  also 
increases. 

3.3.  Dynamic  simulation 

In  this  study,  the  dynamic  simulation  of  three  scenarios  has  been 
performed,  i.e.,  (1)  a  step  decrease  in  acetate  feed  flowrate  from 
2.25  x  10-5  to  1  x  10-5  (m3  h-1 ),  (2)  a  step  decrease  in  acetate  feed 
concentration  from  1.56  to  1  (mol  m-3 ),  and  (3)  periodic  changes  in 
acetate  feed  flowrate  from  2.25  x  10-5  to  1  x  10-5  (m3  h-1 ).  In  each 
scenario,  the  cell  current  density  is  specified  at  2  and  6  (Am-2), 
respectively. 

The  result  of  the  first  scenario  is  shown  in  Fig.  6.  The  acetate 
concentration  is  reduced  with  a  decrease  in  acetate  feed  flowrate 
(Fig.  6a)  as  explained  by  Eq.  (17),  and  it  is  lowered  even  further  if  a 
higher  cell  current  density  is  applied. 

Qc  =  C™  -450/cd,^-  (17) 

The  anode  overpotentials  exhibit  opposite  behaviour  to  the  step 
decrease  in  acetate  feed  flowrate  for  different  levels  of  cell  current 
density  (Fig.  6b).  In  particular,  the  anodic  overpotential  decreases 
when  operating  at  lower  current  density,  and  vice  versa.  Since  the 
cathodic  overpotential  is  not  affected  by  the  change  in  acetate  feed 
flowrate,  the  trends  of  cell  voltage  are  opposite  to  those  of  the  anode 
overpotentials  (Fig.  6c).  As  a  result,  the  power  output  increases 


when  operating  at  lower  current  density  and  reduces  at  higher 
current  density  (Fig.  6d). 

Dynamic  simulation  of  the  second  scenario  is  shown  in  Fig.  7. 
Since  a  decrease  in  feed  concentration  enhances  the  anode  overpo¬ 
tential  (Eq.  (18))  whereas  the  cathode  overpotential  is  not  affected 
by  changes  in  acetate  feed  concentration,  the  cell  voltage  decreases 
as  the  feed  concentration  reduces.  Furthermore,  higher  current 
density  leads  to  an  even  larger  anode  overpotential,  and  thus  a  large 
drop  in  the  cell  voltage  is  observed  (Fig.  7a).  The  resultant  power 
density  is  the  product  of  cell  voltage  and  the  cell  current  density 
(Fig.  7b). 


RT  . 

,?a  =  ^ln 


On  +  VaKdecfx 
k^YacAmfx 


I<Ac 

r  1  (71m  /  Qa) 


(18) 


The  result  of  the  simulation  for  the  third  scenario  is  shown 
in  Fig.  8  with  the  feed  flowrate  switching  from  2.25  x  10-5  to 
lxl0-5m3h-1  at  equal  intervals  of  10  and  2h,  respectively. 
The  resulting  average  power  density  of  the  two  cases  are  1.33 
and  1.31  Wm-2  at  a  cell  current  density  of  1  =  2 Am-2,  and  2.14 
and  2.11  Wm-2  at  7  =  4 Am-2,  respectively.  As  can  be  seen,  at 
7  =  2  Am-2,  the  resulting  average  power  density  of  periodic  feed¬ 
ing  is  1.32  Wm-2,  which  is  higher  than  that  of  the  steady-state, 
i.e.,  1.01  Wm-2.  Similarly,  at  7  =  4 Am-2  periodic  feeding  results  in 
2.11  Wm-2,  i.e.,  higher  than  that  of  the  steady-state,  1.69  Wm-2. 
Furthermore,  the  choice  of  different  intervals  of  flowrate  switch¬ 
ing  does  not  lead  to  a  significantly  different  average  output,  but 
smaller  intervals  tend  to  generate  narrower  variations  which  imply 
a  more  stable  power  output.  The  result  indicates  that  feeding  fuel 
in  a  periodic  fashion  may  achieve  a  higher  average  power  output. 


3.4.  Optimization 

The  optimum  power  output,  subject  to  the  existing  MFC  config¬ 
uration  and  design,  is  studied  in  this  section  by  manipulating  the 
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Fig.  7.  Results  of  dynamical  simulation  of  step  change  in  acetate  feed  concentration  from  1.56  to  1  (mol  m-3):  (a)  cell  voltage;  (b)  power  density. 
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Fig.  8.  Comparison  of  power  outputs  resulting  from  periodic  switches  of  acetate  feed  flowrates  from  2.25  x  10-5  to  1  x  10-5  m3  h-1  with  time  intervals  (a)  At=10h,  (b) 
At =2  h.  Cell  current  density  is  specified  at  J  =  4Arrr2  and  7  =  2  Am-2  for  each  scenario,  respectively. 


Table  4 

Comparison  of  maximum  power  output  at  optimum  operating  conditions  and  existing  condition  for  acetate  MFC. 


flow  rates  of  feed  streams  to  the  anode  and  cathode  compartments. 
In  Scenario  #1,  the  current  density  is  also  allowed  to  vary  so  that  a 
maximum  power  output  is  obtained  at  a  particular  current  density. 
In  Scenario  #2,  optimization  is  performed  to  search  the  maximum 
sum  of  power  outputs  over  a  set  of  equally  spaced  intervals  of  cur¬ 
rent  densities  within  the  range  of  the  experiment.  The  results  are 
presented  in  Table  4  and  Fig.  9,  and  compared  with  the  existing 
experiment.  It  is  found  that  in  Scenario  #1,  the  highest  power  out¬ 
put  2.9  Wm-2  is  achieved  at  current  density  7  =  4.8  Am-2,  but  the 
optimum  operating  condition  only  sustains  proper  operation  up 
to  7  =  6.5  Am-2,  which  is  much  lower  than  the  existing  experiment 
viz.,  7  =  11.5  A  m-2.  In  Scenario  #2,  the  current  density  is  no  longer  a 
manipulated  variable  for  the  optimization,  thus  the  resulting  oper¬ 
ating  condition  (feed  flow  rates)  is  applicable  to  the  entire  range 
of  the  experiment.  The  resultant  maximum  power  output  is  lower 
than  that  of  Scenario  #1,  i.e.,  2.4  W  m-2  at  7  =  7  A  m-2,  but  still  higher 
than  that  of  the  experiment  maximum  of  2.0  Wm-2.  The  result 
implies  that  the  existing  experiment  has  been  operated  at  a  sub¬ 
optimum  condition,  and  reduction  in  the  acetate  feed  flowrate  may 
enhance  the  power  output.  Further  decrease  in  the  flowrate  would 
reduce  the  attainable  current  density. 

3.5.  Simulation  ofGGAfuel  cell 

In  this  work,  the  modelling  method  is  extended  to  a  MFC  fed  by 
artificial  wastewater,  i.e.,  a  solution  of  glucose  and  glutamic  acid 


(GGA).  The  design  configuration  of  the  GGA  MFC  is  similar  to  the 
acetate  MFC  described  in  the  preceding  sections.  The  experimental 
data  reported  in  [27]  is  used  for  comparison  with  the  simulated 
result. 

The  major  difference  between  the  GGA  and  acetate  MFCs  lies  in 
the  bio-electrochemical  reactions  occurring  in  the  anodic  chamber, 


Fig.  9.  Comparison  of  power  densities  among  experimental  measurements  (bold 
line),  by  optimizing  both  feed  flowrate  and  current  density  (dashed  line),  and  by 
solely  optimizing  feed  flowrate  over  entire  range  of  current  densities  (fine  line). 
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Table  5 

Parameters  of  experiment  set-up  and  operation  of  GGA  MFC  model. 
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Symbol 

Description 

Unit 

Value 

T 

Operating  temperature 

K 

306 

Va 

Volume  of  anode  compartment 

m3 

2  x  10-5 

Vc 

Volume  of  cathode  compartment 

m3 

2  x  10“5 

Am 

Area  of  membrane 

m2 

2.4  x  10-3 

dm 

Thickness  of  membrane 

m 

5  x  10“3 

dceH 

Distance  between  anode  and  cathode  in  the  cell 

m 

8.33  x  10-3 

Qa 

Flow  rate  of  fuel  feed  to  anode 

m3  h_1 

2.1  x  10“5 

Qc 

Flow  rate  feeding  to  cathode  compartment 

m3  h-1 

5.7  x  10“4 

CS2 

Concentration  of  dissolved  02  in  the  influent  of  cathode  compartment 

mol  m3 

0.231 

u° 

Cell  open  circuit  potential 

volt 

0.75 

'<2° 

Forward  rate  constant  of  cathode  reaction  at  standard  condition  (from  acetate  MFC  model) 

m12  mol-4  h_1 

3.29  x  10-5 

I<o2 

Half  velocity  rate  constant  for  dissolved  oxygen 

mol  nrr 3 

4.06  x  10“3 

while  the  reduction  reactions  and  mass  balances  in  cathodic  cham¬ 
ber  of  GGA  MFC  are  the  same  as  those  of  the  acetate  MFC.  Since  the 
feeding  fuel  is  a  mixture  of  glucose,  glutamic  acid  and  other  inor¬ 
ganic  compounds,  the  anodic  reactions  are  taken  as  two  oxidation 
reactions  of  glucose  and  glutamic  acid,  namely: 


C6H1206  +  6H20  6C02  +  24H+  +  24e“ 

(RU) 

C5H9N04  +  6H20  -*  5C02  +  NH+4  +  17H+  +  18e“ 

(R12) 

The  two  reaction  rates  (Eqs.  (19)  and  (20))  take  a  similar  form 
as  that  of  acetate  MFC.  k11  and  k12  denote  the  rate  constants  of  the 
corresponding  reactions,  C6  and  C5  are  the  concentrations,  and  I<6 
and  I<5  are  the  half-velocity  rate  constants  of  glucose  and  glutamic 
acid,  respectively.  oq  and  a2  designate  the  charge-transfer  coeffi¬ 
cients  of  the  two  reactions,  respectively.  It  is  noted  that  r\a  (anode 
overpotential)  and  X  (biomass  concentration)  are  the  same  as  those 
used  in  acetate  MFC,  but  their  values  may  be  different  in  the  GGA 

MFC. 

r11=k0exp  (rV^)K6+6C6X 

(19) 

r12=k?2exp(R2r;?Q)K5+5C5X 

(20) 

The  mass  balances  of  five  components  in  the  anode  chamber, 
i.e.,  concentrations  of  glucose,  glutamic  acid,  dissolved  C02,  H+ 
and  biomass,  are  formulated  similarly  to  those  of  the  acetate  MFC, 

namely: 

VaW  =  Qa(C6n-C6)-Anrn 

(21) 

Va^-=Qa(Ci5n-C5)-Amru 

(22) 

Va  =  0a(Ci-nO2  -  Cco2 )  +  6 Amr„  +  5Amru 

(23) 

Va^f=  Qa(9?  -  C„)  +  24 Amrn  +  17 Amru 

(24) 

dX  Qa(X!"-X)  .  v  ,  ,  „  „  v 

17a  ^  ~  j  +AmVAV\/(r  11  +r12)  VaKdAwX 

(25) 

Since  the  cathodic  reaction  in  the  GGA  MFC  is  the  same  as  that  of 
acetate  MFC,  the  rate  expression,  Eq.  (2),  remains  unchanged  in  the 
model  of  the  GGA  MFC.  Similarly,  the  mass  balances  in  the  cathode 
chamber,  Eqs.  (7)-(9),  and  the  relationship  of  the  cell  current  den¬ 
sity  and  flux  of  ions  via  the  membrane,  Eq.  (11 ),  is  used  for  the  GGA 
MFC  model. 

The  charge  balance  of  the  anode  chamber  is  re-formulated  in 
Eq.  (26),  while  the  charge  balance  of  the  cathode  chamber  and  the 
cell  voltage  calculation  remain  unchanged  as  Eqs.  (13)  and  (14), 
respectively. 

Ca%=3600ice„-24Fr1  (26) 


The  experimental  parameters  that  are  different  from  the  acetate 
MFC  are  listed  in  Table  5.  Six  data  sets  of  experiments  [27]  with 
various  combinations  of  artificial  wastewater  flowrates  (Qa)  and 
concentrations  [CAW,  in  the  unit  of  chemical  oxygen  demand,  COD) 
are  used  in  this  study.  The  concentrations  of  glucose  and  glutamic 
acid  in  the  influents  are  in  the  ratio  of  1 : 1  in  all  the  experiments.  The 
range  of  Qa  (ml  min-1 )  is  0.15-0.65  and  CAW  (COD  mg  1_1 )  100-300, 
respectively.  In  particular,  the  six  sets  of  experiments  are:  (a) 
CAW  =  100,  Qa  =  0.35;  (b)  CAW  =  300,  Qa  =  0.15;  (c)  CAW  =  200,  Qa  =  0.65; 
(d)  CAW  =  200,  Qa  =  0.35;  (e)  Qw  =  300,  Qa  =  0.35;  (f)  Qw  =  300, 
Qa  =  0.53.  The  five  model  parameters  (see  Table  6)  are  estimated 
by  using  (a-d)  sets  of  experimental  data  (see  Fig.  lOa-d),  and  the 
model  is  examined  with  (e  and  f)  sets  of  experimental  data  (see 
Fig.  lOe  and  f). 

It  is  shown  that  all  the  simulation  results  generally  reproduce 
the  trends  of  the  experimental  results.  On  the  other  hand,  relatively 
large  deviations  are  observed  when  the  wastewater  concentration 
(CAW)  or  flowrate  (Qa)  is  at  the  lower  ends  of  the  respective  ranges  in 
Fig.  10a  (CAW  =  100  COD  mg  H )  or  Fig.  10b  (Qa  =  0.15  ml  min-1 ).  The 
larger  discrepancies  may  be  caused  by  inconsistent  experimental 
results  [27],  which  showed  that  the  I-V  curves  obtained  from  the 
two  operation  conditions  significantly  deviate  from  the  others  oper¬ 
ated  at  higher  CAW  or  Qa.  The  discrepancies  may  also  arise  from  the 
fact  that  the  model  is  unable  to  deal  with  an  GGA  MFC  operating  at 
extreme  conditions,  thus  further  improvement  is  required. 


Table  6 

Estimated  parameters  of  GGA  MFC. 


Symbol 

Description 

Unit 

Resultant  values 

fc?i 

Forward  rate  constant  of  glucose  reaction  at  standard  condition  (maximum  specific  growth  rate) 

mol  m  2  h_1 

5.94  x  lCH 

I<6 

Half  velocity  rate  constant  for  glucose 

mol  m3 

3.68  x  10-5 

k°u 

Forward  rate  constant  of  glutamic  acid  reaction  at  standard  condition  (maximum  specific  growth  rate) 

mol  nrr  2  h_1 

7.35  x  10-3 

K 5 

Half  velocity  rate  constant  for  glutamic  acid 

mol  nrr3 

1.42  x  lCH 

ai 

Charge  transfer  coefficient  of  glucose  oxidation  at  anode 

Dimensionless 

0.760 

ol2 

Charge  transfer  coefficient  of  glutamic  acid  oxidation  at  anode 

Dimensionless 

0.753 

fi 

Charge  transfer  coefficient  of  cathode 

Dimensionless 

0.740 
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Current  Density  (A  m'2) 


Current  Density  (A  m'2) 

Fig.  10.  Comparison  of  experimental  (diamonds)  and  simulated  (triangles) 


Current  Density  (A  rrf2) 


Current  Density  (A  nrf2) 

results  for  a  GGA  MFC,  (a-d)  fitting  experiments;  (e  and  f)  model  validation. 


4.  Conclusions 

The  present  work  describes  a  method  for  modelling  the  micro¬ 
bial  fuel  cell  (MFC).  The  basis  of  the  method  is  mass  and  charge 
balances  with  the  integration  of  bio-electrochemical  reactions.  It 
is  worthwhile  to  note  that  the  model  describes  the  experimental 
findings,  i.e.,  virtually  no  protons  are  transported  in  the  MFC  and 
electro-neutrality  is  mainly  sustained  by  the  transport  of  cations 
instead  of  protons.  The  results  of  this  modelling  explain  the  trends 
in  the  experimental  data  under  a  steady  or  quasi-steady  state,  and 
give  insights  into  how  various  parameters  affect  the  power  output. 
In  particular,  the  cathodic  reaction  may  be  the  most  significant  fac¬ 
tor  in  limiting  the  performance  of  MFCs,  which  is  in  agreement  with 
other  published  studies.  Steady-state  simulation  demonstrates  that 
a  decrease  in  the  fuel  feed  flowrate  may  lead  to  a  boost  in  power 


output,  but  at  the  expense  of  reducing  the  attainable  current  den¬ 
sity.  Furthermore,  the  simulation  also  shows  that  an  increase  in 
fuel  concentration  is  favorable  in  terms  of  power  output  and  also 
allows  the  MFC  to  operate  at  a  higher  current  density.  The  results 
of  dynamic  simulation  suggest  that  a  periodic  change  in  fuel  feed 
flowrate  can  yield  higher  average  power  output  compared  with 
steady  state  operation,  which  provides  guidance  for  more  efficient 
MFC  designs  and  operational  strategies. 

While  the  model  is  developed  based  on  an  acetate  MFC  with 
a  typical  two-chamber  design,  the  present  study  shows  that  the 
method  can  also  be  applied  to  a  MFC  fed  with  other  fuels  such 
as  a  solution  of  glucose  and  glutamic  acid.  It  is  noted  that  since 
the  lump  parameter  method  is  used  to  formulate  the  model,  some 
of  the  details,  such  as  those  occurring  in  the  biofilm,  cannot  be 
predicted  from  the  model.  Nevertheless,  models  of  MFCs  with 
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various  designs  or  operating  conditions  can  be  readily  developed 
using  this  relatively  simple  method,  which  serves  as  a  starting 
point  for  the  development  of  more  comprehensive  MFC  models 
that  address  more  complicated  issues  relevant  to  MFC  perfor¬ 
mance. 
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